Relative Critical Points 



(N 



00 



Debra LEWIS ^ 

^ Mathematics Department, University of California, Santa Cruz, Santa Cruz, CA 95064, USA 
E-mail: lewis@ucsc.edu, 

Abstract. Relative equilibria of Lagrangian and Hamiltonian systems with symmetry are 
critical points of appropriate scalar functions parametrized by the Lie algebra (or its dual) 

■ of the symmetry group. Setting aside the structures — symplectic, Poisson, or variational — 
. generating dynamical systems from such functions highlights the common features of their 

construction and analysis, and supports the construction of analogous functions in non- 
, Hamiltonian settings. If the symmetry group is nonabelian, the functions are invariant 

^w' ■ only with respect to the isotropy subgroup of the given parameter value. Replacing the 

parametrized family of functions with a single function on the product manifold and ex- 
tending the action using the (co)adjoint action on the algebra or its dual yields a fully 
invariant function. An invariant map can be used to reverse the usual perspective: rather 
than selecting a parametrized family of functions and finding their critical points, condi- 
tions under which functions will be critical on specific orbits, typically distinguished by 

■ isotropy class, can be derived. This strategy is illustrated using several well-known mechan- 
, ical systems — the Lagrange top, the double spherical pendulum, the free rigid body, and the 

Riemann ellipsoids — and generalizations of these systems. 
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O ; 1 Introduction 

d. 

The analysis of steady motions of dynamical systems with symmetries has played a pivotal role 

■ in the development of geometric mechanics, and the identification of such motions is central to 
the design and interpretation of mathematical models arising in a broad range of disciplines. 
Relative equilibria of Hamiltonian or Lagrangian systems with symmetry can be characterized 

^ I as critical points of functions, often parametrized by elements of the algebra (or its dual) of the 

^ ■ symmetry group. Steady motions with nontrivial isotropy frequently serve as the starting point 

for both conceptual and detailed analyses of systems with symmetry. (See, e.g., [MllSl].) Highly 
symmetric states are often more intuitively accessible than more general motions; once some 
families of steady motions are known, a variety of techniques can be applied to locate additional 
trajectories with more subtle structure. (See, e.g., 134 1 [T7 1 122j). Many classical and modern 
applications, e.g., Lagrange tops and isotropic hyperelasticity, involve both spatial and body 
symmetries; for such systems, physically symmetric steady motions such as sleeping tops and 
steadily rotating spheroidal configurations, are readily predicted and visualized. The relative 
equilibria for such systems frequently include families of quasi-periodic motions bifurcating from 
families with continuous isotropy subgroups. The evolution of the model now known as the 
Riemann ellipsoids is a striking example of this process (see f|5l |38t El [12] ) . 

Traditionally, selection of an appropriate dynamical model is the first step in the mathemat- 
ical analysis of a system; after the model is specified, its relative equilibria are identified and 
classified. When constructing a model of data for which the key governing phenomena are not 
known a priori, but for which certain symmetries are believed to be crucial or certain steady 
motions have been observed, we may seek dynamical systems capturing these features. In this 
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situation, it may be desirable to identify all systems within a diverse family that have the rele- 
vant features. Characterizing such systems as elements of the kernel of an appropriate operator 
yields computationally tractable conditions that can be combined with standard data-fitting 
techniques. 

At present, there appears to be no overall best method for the analysis of relative equilibria — 
each approach performs well for important classes of problems, but is inapplicable or awkward 
to implement in many equally important situations. Riemann's analysis of the existence and 
stability of steady affine motions of a very simple liquid drop model employs techniques that are 
recognizable to the modern reader as precursors of several of the methods described below |38pi2j 
(see different strategies are used for each of the classes of steady motions, exploiting the 
special features of each class. Thus it seems advantageous to emphasize the common features of 
this family of techniques, avoiding wherever possible the invocation of specialized structures that 
are not essential to the steps under consideration. We will focus strictly on the critical points 
of the functions related to group actions on manifolds, ignoring any underlying variational, 
symplectic, or Poisson structure. Our intent is both to organize a versatile tool kit and to 
distinguish components of traditional relative equilibrium analyses that make explicit use of 
these geometric structures from those that utilize only the group action. 

Relative equilibrium methods that make use of an equivariant momentum map typically either 
explicitly restrict attention to a single momentum level set in some stages of the analysis or use 
a Lagrange multiplier approach to implement a constrained criticality test. If the symmetry 
group is nonabelian, both of these approaches break the symmetry of the original system unless 
the momentum value under consideration has full isotropy, i.e. the isotropy subgroup with 
respect to the coadjoint action is the entire group. The functions resulting from such approaches 
are invariant only with respect to the isotropy subgroup — in some situations, e.g. symplectic 
reduction, the full group does not act on the new manifold of interest; in others, e.g. the energy- 
momentum method, the original manifold is retained, but the original function is replaced by 
one that is no longer invariant with respect to the full action. In the latter situation, 'rigid' 
variations tangent to the group orbit play a central role in the analysis, but the information 
captured by such variations cannot be accessed by applying invariant function theory to the 
new reduced-symmetry function. 

To maintain the full group invariance, we regard the (dual) algebra element as a variable, 
not a parameter, and extend the group action to the product manifold using the (co)adjoint 
action to obtain a fully invariant function. This allows the use of invariant function theory for 
the full symmetry group and facilitates systematic organization of point-algebra element pairs 
by isotropy class. We only require criticality with respect to the original variables. However, 
utilizing the invariance of the function, we can replace variations tangent to the group orbit in 
the manifold with variations tangent to the (co) adjoint orbit in the algebra, obtaining the so- 
called rigid criticality conditions. This strategy can reveal purely symmetry-dependent features 
of the system. 

Steady affine motions of a self-gravitating fluid mass were studied by some of the greatest 
mathematicians of the 18th and 19th centuries, including Newton, Maclaurin, Jacobi, Dirichlet 
|llj . Dedekind [10], and Riemann [38]; the techniques introduced in the analysis of this system 
foreshadow many key components of geometric mechanics, bifurcation theory, and Lie group 
theory. We briefly explore the relative equilibrium conditions for some classes of configuration- 
generator pairs identified by Riemann |38j . 
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2 Relative equilibria as critical points 

The general technique for variational characterizations of relative equilibria can be briefly de- 
scribed as follows: Consider a manifold P with a G action and a G-equivariant vector field X. 
Given a point pe £ P and an element ^ of the Lie algebra q of G, we wish to determine conditions 
on pe and ^ under which the solution of the initial value problem p = X{p) and p(0) = pe is 

=exp(t^)-pe. (1) 

Such trajectories are known as relative equilibria or steady motions. Differentiating ([T]) with 
respect to time, we see that a point p is a relative equilibrium with generator ^ of a dynamical 
systems with vector field X iff X{p) = $,p{p), where the vector field 

^pip):=iexp{t^)-p\t=o (2) 

is the infinitesimal generator associated to the algebra element ^. At a relative equilibrium p^ 
with continuous isotropy subgroup the generator $, is not unique; $, may be replaced by ^ + ^ for 
any element C of the isotropy algebra of pe- Both continuous and discrete isotropy play a crucial 
role in the analysis of relative equilibria. (See \23\ I24j for detailed discussions of this issue.) 



2.1 Phase space critical point characterizations 

We first briefly review three well-known characterizations of relative equilibria as critical points 
of functions on the full phase space, beginning with the most traditional and explicit, and moving 
in the direction of increasing generality. 

• Energy-momentum method (Lagrangian version) 

Consider a Lagrangian system with configuration manifold Q, symmetry group G acting 
on Q, and Lagrangian L : TQ M invariant under the lifted action of G on TQ. The 
energy of the system is given by the function E : TQ M defined by 

Eiv,) := FLK)K) - L{v,) = £ L((l + e) v,) - e Liv,)^^, . (3) 

Here ¥L : TQ T*Q denotes the Legendre transformation 

Noether's Theorem gives the momentum map J : TQ q* satisfying 

j^{v,)=¥L{v,){^Q{q)). 

The vector field determined by the Euler-Lagrange equations with Lagrangian is the 
infinitesimal generator ^^"q © associated to ^ € g and the lifted action of G on TQ. Thus 
Vq is a relative equilibrium with generator ^ iff Vq is a critical point of E — j^. (See, e.g. 

mm-) 

• Energy-momentum method (symplectic version) 

Let {P, oj) be a symplectic manifold and G be a Lie group acting on P by symplectomor- 
phisms. If there is a momentum map J : P q* for the G action, i.e. a map such that 
for any ^ € g the scalar function given by j^{p) := J{p) ■ has Hamiltonian vector field 
equal to the infinitesimal generator of ^. p is a relative equilibrium for a Hamiltonian 
iff p is a critical point of H - j^. (See, e.g., [D [H [25l [261 [271 [28].) 
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• Energy-momentum-Casimir method 

Let P be a Poisson manifold with Poisson bracket { , }, and let G act on P by Poisson 
maps. A non-constant smooth function C : P — > M is said to be a Casimir if {F, C} = 
for all smooth functions F on P. If has momentum map J : P — ?• g*, i.e. {F, J^} = i,p[F] 
for all smooth functions F on P, then p is a relative equilibrium with generator ^ iff p is 
a critical point H — + C for some Casimir C. (See, e.g., [29l [30l [32| [18].) 

Many Poisson manifolds and noncanonical symplectic manifolds in geometric mechanics arise 
through the process of reduction. In symplectic reduction, the reduced phase space associated 
to the momentum value /x S g* is defined as the quotient of the fi level set of the momentum 
map by the subgroup of group elements preserving /x. The symplectically reduced phase space 
has an induced symplectic structure; the reduced equations of motion are defined using this 
symplectic structure and the induced Hamiltonian. (See, e.g., [33l[32l[Tll36].) Poisson reduction 
proceeds as follows: Given a Poisson manifold P with a free and proper G action by Poisson 
maps, the Poisson bracket on the quotient manifold M := P/G with projection map vr : P — )• M 
is determined as follows: for any smooth real-valued functions / and h on M, require that 

{f: hj^oTT = {f oTr,hoTr}p. 

(See, e.g., [301 [19].) Complications arise in reduction if the group action is not free. Continuous 
isotropy leads to conical singularities of the momentum level sets [4], and singularities of the 
quotient manifold arise even in the presence of discrete isotropy. Various approaches to ad- 
dressing these difficulties have been developed: Singular reduction has been implemented using 
stratified symplectic spaces [H [33] and spaces of invariant functions [21 [3] . In situations where 
the isotropy is due to the interaction of multiple components of a product group action — e.g., 
body and spatial symmetries of a rigid body or continuum system — one or more components 
of the group can be excluded from the reduction process and subsequently incorporated in the 
analysis via application of the energy- momentum method to the reduced manifold. 

Reduction methods are valuable analytic tools in the analysis of relative equilibria, but are not 
always convenient to implement numerically or symbolically. On a fully symplectically reduced 
phase space, relative equilibria are critical points of the Hamiltonian; hence there is no need to 
match appropriate Casimir functions to steady motions. However, there no longer is a single 
phase space for the entire system. Rather, there is a distinct reduced phase space for each value 
of the momentum and these spaces can change dramatically as the momentum value is varied. 
The analysis of these changes can be quite complicated, and require a substantial amount of 
topological and geometric sophistication, even for very low dimensional, simple systems. In the 
presence of nontrivial isotropy, these complications are even greater. The reduced phase space 
is no longer a manifold, but instead is a stratified symplectic space. The variational analysis of 
relative equilibria is carried out stratum by stratum. If partial reduction is performed, resulting 
in a noncanonical structure with a momentum map associated to the remaining group action, 
the inequitable treatment of the components of the symmetry group can obscure the relationship 
between generators and steady motions, particularly with points which had nontrivial isotropy 
in the original manifold. 

2.2 Relative equilibrium characterizations on the configuration manifold 

Lagrangian systems and their canonical Hamiltonian analogs on (co)tangent bundles with lifted 
actions have alternative critical point characterizations. The key observation motivating the 
reduced energy momentum (REM) method and related approaches is very simple: due to the 
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structure of the infinitesimal generators determined by the lifted action, the velocity of a relative 
equilibrium Vq £ TQ with generator ^ satisfies Vq = ^giq)', the associated canonical Hamiltonian 
systems on T*Q has relative equilibrium of the form ¥L{^Q(q)), where ¥L denote the fiber 
derivative of the Lagrangian L. Thus once an algebra element ^ has been selected, there is no 
need to explicitly solve for the velocity — all relative equilibria with generator ^ are uniquely 
determined by the base point q. Given an element ^ of the Lie algebra q (respectively element 
/X of g*), we can use a variational criterion to derive conditions on configurations guaranteeing 
that they are basepoints of relative equilibria with generator ^ (respectively, determine a state 
with momentum map value fi). 

Smale's analysis of relative equilibria of simple mechanical systems utilizes a characterization 
of points in the cotangent bundle generated by infinitesimal group motions as minimizers of 
kinetic energy within a given level set of the momentum map [441 [45] . This approach yields a 
variational characterization of relative equilibria of such mechanical systems on the configuration 
manifold Q. Assume that a simple mechanical system with configuration manifold Q has energy 
E{vq) = ^ \\vq\\^ + V{q). Given /2 S g*, the amended potential : Q M associated to n is 
given by 

V^{q) = \^l■l{qr'^l-V{q), (4) 

where \{q) : g — )■ g* is given by 

{l{q)i)-V-= {^Q{q),VQ{q))q 

and fi = I((?)^. (Since the action on a simple mechanical system is free, l{q) is invertible.) This 
explicit construction is linked to Smale's constrained energy minimization characterization as 
follows: if we define : Q ^ J^^(/x) by 

^^/x(g) := (I(9)"V)q(9), then Vfj, = EoVfj,. 

Critical points g of correspond to relative equilibria v^{q); local minima (modulo symmetries 
of the equilibrium momentum /i) of correspond to (formally) stable relative equilibria |28[ W\\ . 

The strategy underlying the REM method for simple mechanical systems can be extended to 
more general functions and symmetries, relaxing the assumptions that the Lagrangian have the 
form 'kinetic — potential' and that the action be free. Given q € Q and ^ G g, ^q(q) is a relative 
equilibrium of the Euler-Lagrange system with G-invariant regular Lagrangian L : TQ — > M 
iff g is a critical point of the locked Lagrangian := L o ^q] see |23[ I24j . One of the key 
advantages of the (relative) energy-momentum method relative to singular symplectic reduction 
techniques is the ease with which continuous isotropy is managed. The locked Lagrangians have 
no singularities, even at points with jumps in isotropy; since the momentum level sets do not 
play an explicit role in the analysis, possible singularities of the level sets are not a concern. 

3 Critical points of invariant functions 

Motivated by the considerations touched on in the previous section, we consider critical points of 
families of functions parametrized by a Lie algebra or its dual, ignoring any underlying geometric 
structure beyond a group action. If the symmetry group is abelian, we can use parametrized 
families of real-valued group invariant functions. If the group action is nonabelian, we replace 
the parametrized family of functions with a single function depending on both the given manifold 
and the algebra (or its dual) of the symmetry group, requiring that this function be invariant 
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with respect to the action on the product manifold determined by the original action and the 
(co)adjoint action on the algebra or its dual. We focus on situations in which these functions 
can be expressed as the composition of an invariant map to an intermediate manifold with some 
scalar function on that manifold. This facilitates the identification of common properties of a 
wide class of functions sharing common symmetry properties. 
We begin with two familiar, elementary examples. 



3.1 The Lagrange top 

The heavy Lagrange top is one of the simplest mechanical systems with a nonlinear configuration 
manifold and a nonfree group action. Hence we take this system as our first example of the 
strategies described above. The total energy of the heavy top is invariant under spatial rotations 
about the axis of gravity and body rotations about the axis of symmetry of the top. See, e.g., 
[391 [231 [21] • 

A heavy sleeping top, i.e. a top rotating at constant speed under the influence of gravity, 
with the spatial axis of symmetry aligned with the direction of the force of gravity, is fixed by 
the combination of a spatial rotation through an angle 9 and a body rotation through the angle 
—9. Note that spin and precession are not uniquely determined for a sleeping top — spin about 
the axis of symmetry and precession about the axis of gravity are indistinguishable when those 
axes are aligned. This continuous isotropy significantly complicates many geometric approaches. 
Here we briefly present an investigation of the relative equilibria of the Lagrange top using the 
locked Lagrangian; see [23] . 

The Lagrange top has Lagrangian 

L{A^) = ^n'^m- gmejAa, (5) 

where 

• I = /i 1 + (/a — /i)<T<T^ is the inertia tensor of the top, for positive constants Ii and I3, 
and unit vector cr determining the axis of symmetry of the top 

• 5^ 63 is the force of gravity 

• mcr is the center of mass. 

The Lagrangian ([5]) is invariant under the action of x on SO{3) by 

{9,(p) ■ A = ex-p{9 es)Aexp{—(pa), (6) 

corresponding to spatial rotations about the vertical axis 63 and body rotations about the axis 
of symmetry a. 

The locked Lagrangian L(^^^^^-) : S0{3) M determined by L and the algebra element 
(C^) Cr) S ~ is a quadratic polynomial in the invariant l : 50(3) — )• [—1, 1] given by 

i{A) := eU^, (7) 

with differential 

dL{A){Ar]) = eJrjAcr = {e-s x Aa)'^r], 

where ^ : ^ {3 x 3 skew-symmetric matrices} is given by ^ = x for all ^ € M'^. Thus the 
critical points of i are the 'sleeping' states, with Aa = ±63; the critical point set is 1, 1}), 
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the preimage of the endpoints of the codomain of l. Sleeping states are fixed by the S action. 
All other points have trivial isotropy. 

The locked Lagrangian satisfies L(^^^g^) = ^(5^.5^) o where 

ki^&)i^) ■■= hfiih - h)x^ - i^e &h+9 m)x + \{e,h + Cr^s)- 

Clearly, sleeping tops are critical points of L^^^^^^) for any parameter values. 'Tilted' states, i.e. 
elements of 1)), are critical points of L(g^ ,f^) iff 

e,{h-h)i{A) = i,irh + 9m. (8) 

For any A E 50(3) and nonzero value of there is a unique value of such that A is a critical 
point of and hence determines a relative equilibrium with generator (^£,^r)- 

The only true equilibria are the sleeping states. Note that £(o,o)(^) = —gmx is linear; hence 
the critical points of L[ofi) are precisely those of l. The Lagrangian L is invariant with respect 
to negation of the velocity, corresponding to the time reversal symmetry of the dynamics; this 
invariance induces the 'symmetry' that two algebra elements differing only by sign determine the 
same locked Lagrangian. The origin is the unique fixed point of the Z2 action on M^, ^ 1— )■ 
Of course, in the case at hand, it is hardly surprising that the conditions for true equilibrium are 
significantly simpler than those for relative equilibrium; however, we shall see that in systems 
with more elaborate symmetry groups, consideration of the relevant group action on the algebra 
and classification of the algebra elements with regard to isotropy can be a powerful tool in the 
analysis of relative equilibria. 

Remark. An upright state, with l{A) = 1, has isotropy algebra span{(l, 1)}; setting (^^,^r) = 
(A, A — ^) in ([8]) sheds light on the relationship between the upright relative equilibrium and 
nearby tilted states with trivial isotropy. Solving for A yields 

_ M ± ^ei'i - ^gmjh + (ii - /3)^(^ 

2(/3 + ih - l3)i{A)) 

Thus, setting i{A), we see that the family of solutions ([8]) intersects the family of sleeping tops 
throughout the range of angular velocities (,^1^ > Agmli for which the sleeping tops are stable. 
Corresponding to every speed for which the sleeping top rotates stably, there is a family of 
tilted, precessing motions whose combined spin and precession equal the angular velocity of the 
sleeping top [H]. Some approaches, particularly those involving stratified spaces or reduction 
with respect to a subgroup, can obscure this simple but significant relationship between the 
generators of sleeping and tilted states. □ 

The division of critical points into sleeping and tilted states holds for any parametrized 
function of the form F = f o l for some function / : M ^ R. Sleeping states are critical points 
for any function /; additional critical points of F are states A such that l{A) is a critical point 
of/. 

3.2 The double spherical pendulum. 

Our second mechanical example with an abelian symmetry group is the double spherical pen- 
dulum. In this example, discrete isotropy subgroups play a central role in the classification of 
the relative equilibria. 
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Let Qj € denote the position of the j-th pendulum, with mass rrij and length ij. Let 
g denote the strength of gravity and ^ G M denote the generator, i.e. the (vertical) angular 
velocity of the pendula. The Lagrangian for the double spherical pendulum is 

HQi,(12) = ^ Mlf + ^ IIQi +q2f -963 ((mi +7712)91 +771292)- 

The relative equilibria of the double spherical pendulum are studied using a combination of 
Lagrangian reduction and blow-up techniques in [21]. Here we again use the locked Lagrangian 

L(iqi,q2) = T ("^1 11^3 X Qif +7772 ||e3 X (Qi +92)11^) " 5 63 (("il + "^2)91 + "72^2) • (9) 



The locked Lagrangian ([9]) is invariant under the diagonal action of 0(2) consisting of rota- 
tions about the vertical axis and reflections across vertical planes. We first formulate critical 
point conditions for a broad class of functions on S'^ x S"^ invariant under the action de- 
scribed above, then specialize those conditions to the locked Lagrangian for the double spherical 
pendulum. We define the invariant map t : x — t- given by 

L{q) := (ea • Qi, eg • 92, 9i • IP392), (10) 

where P3 denotes projection into the horizontal plane, and consider invariant functions of the 
form F = f o L for some function / : — )■ M. Note that if {91,921^3} linearly dependent, 
then 

4 = {I -4) {I -4)- (11) 

The critical point behavior is determined by the isotropy of 9, which is in turn determined 
by the dimension of span{9i, 92, 63}: 

Proposition 1. 

• // 9i and 92 are both vertical, then q has full isotropy and q is a relative critical point of 
any invariant function. 

• ^/ 92; ^3} spans a plane, then 9 is fixed by reflection across that plane. 

// = 1 7^ \Lji\, where {j,j'} is a permutation 0/ {1, 2}, 9 is a critical point iff 

df , , Of , , 

Oij' Ois 

If 7^ 1 7^ I'-j'l, 9 is a critical point iff 
df Of 

~ ^"^^^ dT-^^'' ^ '^^''^'di^^'^^' -^"^ J = l;2 and i = i{q). (12) 

• // {9]^, 92, 63} is linearly independent, then 9 has trivial isotropy and is a critical point iff 
L{q) is a critical point of f. 

Proof. 

If we define pj := x e^, j = 1, 2, then 
Pi 

di{q) =1 P2 I . (13) 

9i X IP392 Q2 X F391 
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If Qj is vertical, then = 1, = 0, and pj = 0. If both qi and 92 ^-^^ vertical, then 
di{q) = 0; hence q is a critical point of any invariant function. If q^^ is vertical, but is not, 
then 



di[q) 



Since p2 G ^Qj*^^ ^o'^ J = 1)2 in this case, range[(i/-(q)] = span{e2,e3}. The case in which 
|i2| = 1 / kil is analogous. 

If neither pendulum is vertical, we can take 

Pi W W Pi X Qi W 

y ' V P2 y ' I y ' V P2 X 92 

as a basis for the tangent space. With respect to this basis, di{q) has matrix representation 

IPif \ 
Ilp2f ^ ' where r = 63 (q^ x 92)- (14) 

-Lii-i -L2L3 T -T J 

Hence dt^q) has full rank iff r 7^ 0, i.e. {Qi,Q2i^3} is linearly independent. (jl2p now follows 
directly from consideration of the nullspace of the transpose of ([140 . □ 

We now use Proposition [1] to find the relative equilibria of the double pendulum. As in the 
previous example, the locked Lagrangian ([9]) for the trivial generator is linear in the invariant: 

Loiq) = -c/("ii^i,m2^2,0)^t(q). 

The only critical points of Lq are the critical points of l, namely the vertical pendula. To 
streamline the expression for the locked Lagrangians with nonzero generators, we introduce the 
parameter ratios 

. h rn2i g 
^■=T' ^ •= \ ' ^^'^ ^ •= 72- 

£1 mi + 771-2 ? 

For ^ 7^ 0, the locked Lagrangian satisfies = ^^£i(mi + 1712)/ o l for 

f{i) := ^1 - il + mi (1 - lI) + 71113) --f{ii+mi2), (15) 
The function (llSp has differential 

1 - ti - 7 

df{i)= I m(^(l-.2)-7) 
m 

Thus, since ^ = m 7^ 0, Proposition [1] implies that if q is a critical point of the locked 
Lagrangian, then q^ and ^2 in a common vertical plane; either both pendula are vertical or 
neither is vertical. In the case that neither is vertical, we can introduce the expression 

X := — ^-2, which satisfies = 1 (16) 
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when {91,^21^3} is dependent. The criticality conditions (|12p then take the form 

m(£(l - 12) - 7) = ^('') = ''2^^^'') = ''T'^a^^ j 

We can decompose the analysis into subcases determined by the relationships between the 
parameters i and 7: 

.7^1/1 

The criticality conditions PT|) have the unique solution 

1-7 £-7 
^1 = and .2 = 

Substituting these values into the righthand expression in (116p determines a sixth order 
polynomial in x. whose roots determine candidate values of i\ and 1% only those values 
for which < 1, j = 1,2 are admissible. Note that this polynomial can be regarded as 
a quadratic polynomial in 7, with k treated as a parameter if it is preferable to solve for 
the angular velocity ^ as a function of i\ and l^- 

• 7 = 1. 

Either 

,2 

ti = and (£(1 - L2) - \f = 



4 



or 



7 , 



-X = ?n ^ 1^2 = 7i > and 

I — m \m 



Either 



or 



2 ("^'-i)^ 

i2 = and (1 — ii - 7) = 2" 



-x = ^/m, ^1 = X' a-'^d (^^2)^ = 1 - 



3.3 Adequate invariant maps 

Analysing classes of invariant functions, rather than a single function or parametrized family of 
functions, can shed light on the role of symmetry in determining the critical point conditions(See, 
e.g. [13l[16l[35l[36l[371[8].) This can be particularly valuable in situations where the appropriate 
parameterized models are not known a priori and choices must be made on the basis of a 
range of desired behaviors. We outline a possible approach that can be used to symbolically 
or numerically search for functions having specified critical point and isotropy relationships. 
Generalizing the approach of the previous section, we consider invariant functions constructed 
as the composition of an invariant map and scalar functions on the codomain of that map. If 
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t : -P — > Q is G-invariant and F = f o i for some / : Q — > M, then p is a critical point of F iff 
(^Li'p)/ S ker [d* t] . By characterizing these kernels, we can obtain critical point conditions for 
such invariant functions. 

Generating sets of invariants are known for many important group actions, and can be used to 
construct an invariant map t : P Q. (See, e.g., 120 1 [21].) If ^ generating set is not known 
a priori and cannot be readily determined, more limited criteria for diversity of an invariant 
map L : P ^ Q can be used. Such a map may not capture the full range of possible invariant 
function behavior, but can still be used to study very general classes of invariant functions. We 
describe a possible criterion — that the rank of the differential of the map be as large as possible 
at all points — that is relatively easy to verify in many cases, and then illustrate implementation 
of this test in several examples. 

Let : P P denote the action of € G on a manifold P, i.e. ^g{p) = g • P- If -F : P — )■ M 
is G-invariant, then for any ^ G 0, 

^p[^] = l^°'^'exp{U)k=o = f = 0; 

hence for any p €z P, range [dpF] C (g • p)^, the annihilator of the tangent space to the group 
orbit through p. There are additional restrictions on the differentials of invariant functions at 
points with nontrivial isotropy. If 5 € Gp, the isotropy subgroup of p, then 

dF{p) O dp<l>g = dF{^g{p)) O dp^g = d{F O ^ g) {p) = dF (p) 

implies 

range [dp^g - 1] C kei:[dF{p)], i.e. dF{p) G ker [{dp^g)* - 1] . 

Applying this argument to F = foi for a G-invariant map i : P Q and arbitrary differentiable 
real-valued function / on Q, we see that 

range [d*pL] C (g • p)^ n {CigeGp ker [d*p<^g - l]) = (g • :p U (u^gCprange [dp<^g - 1]))^ . 

If 

rank [d^i] + dim [g ■ pU (u^gGprange [dp^g — 1])) = dimP (18) 

for all p G P, we will say that the map l is adequate. 

The following examples illustrate the process of verifying adequacy in some simple and fa- 
miliar settings. 

Examples. 

• The usual action of 0(3) on W^. 

i : M'^ — > M given by t(x) := ||x||^ is adequate, since d^L is rank one and g ■ x is two 
dimensional if x 7^ 0, while range [do*J'-i ~ 1] = range [—2 1] = implies 

UAGGo=0{3)range [do^A - 1] = M^. 

• acting on 50(3) hj 9 ■ A = RgA, where Rg implements rotation through 9 about a 
fixed axis cr ^ S^. 

The invariant map l : 50(3) — t- 5^ given by i{A) := A'^cr, with differential 

dL{A){Ar]) = -rjA^cr = t{A) x 77, 

is adequate, since range [dL{A)] = r^(^)52 for all A G 50(3). This example illustrates the 
possible advantages of regarding i as a map into a manifold, rather than a collection of 
invariant scalar functions. 
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• xS^ acting on SO (3) by 

In H'6.1\ we saw that ^ is a critical point of the invariant map ([7]) iff ^ is a sleeping state, 
with Ao- = ±63 and Ga = {e,e):ee S^}. 

dA^^eM^"^) = Ree-iVAR-ea- = ReeaVReAo-A = RoeaVReAezA = ReesijA 
imphes that 

range [^^^(e,^) — l] = spanjeiA, 62^} if d^l-rrk, A; G Z. 

Combining this with q - A = span {63^,^0^} = span {637!} if Act = ±63, we see that in 
this case 

(0 • ^) U (Ueg^irange - l]) = ^4^0(3) 

at a sleeping state. 

If A is not a critical point of l, then di{A) is rank one and A has trivial isotropy. It follows 
that dimg • A = dim (S^ x S*^) . Hence rank[(ii(yl)] + dimg • j4 = 3 at these points as well, 
and i is adequate. 

• 0(2) acting on S"^ x 5^ as in gS^l 

If we regard Tg.S'^ as a subspace of M^, then 

where P3 denotes orthogonal projection onto the 61-62 plane. If both and ^2 
vertical, then range [dg<^^ -l]=Tg (52 X 52). If and 93 lie in the same vertical plane, 
then Proposition [1] implies that rank[(i*i] = 2; since 

range [dq<^^ - 1] 5 span{(63 x q^, 0), (0, 63 x ^2)} 

in this case, (jlSp is satisfied at these points. If spanlg^, ^2) ^3} = M.^, then rank[(i*t] = 3 
and • q is one dimensional. Hence (jlSp is satisfied at all points and i is adequate. 

• 0(3) x 0(3) acting on GL(3) by {U, V) ■ A = UAV^. 

l{A) := (trace C, \ ((trace Cf - trace C^) , det C) , where C := A. 

Equivalently, setting io = 1) 
3 

det(0- Al) = J]]ij(A)(-A)3-^' VAgM. 

i=o 

We have 

rank[dyi;.] = # distinct singular values of A. 

If A has three distinct singular values, then the isotropy subgroup of A is discrete, 3 • A is 
six dimensional, and the differentials of the components of i are linearly independent. 

If A has only two distinct singular values, then there are matrices U,V G 0(3) and scalars 
Q, d such that 

A = U {al + aesej) V'^ . 
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If we set u := f/e3 and v := det{UV) Ve^, then the isotropy subgroup Ga of A is 

GA = {{Reu,Re.):0 gS^] , 

where denotes rotation about x through the angle |x| for any nonzero x G M'^. If we 
define ipe ■ M^x^ ^ rSxs 

MB) -ReuBRj^-B, 

then [g- A)U (Ugg^irange [V'e]) has orthogonal complement spa.n{UV^ , uv'^} with respect 
to the inner product {B,B) = tiace{B'^ B). If T^GL{3) is identified with W^^^ by means 
of this inner product, then 

range [d^t] = span{UV^ , uv^} 

as well. Hence (llSp holds for A with multiple singular values. 
If A has exactly one singular value, then A € 0(3) and 

Ga = {{AU,UA) : U € 0(3)}. 
The linearization of the action of an element of the isotropy subgroup satisfies 

dA^(AU,UA)(^^) = AU{VA){UAf = AUV{AUfA = (AdAuV) A. 
Hence if pA denotes right multiplication by A, then 

dAPA^ ° {dA^g - 1taGL(3)) ° dlPA = AdAf/ - 1r3x3 

and thus 

dApJ^ (UgeG^range [dA<^g - 1taGL(3)]) = U^go(3)range [AdAu - Irsxs] 

= {traceless 3x3 real matrices}. 

If j4 € 0(3), then range [d'^i] = spanjtrace o dAP~^^}- Hence (jlSp is satisfied for A with a 
unique singular value and l is adequate. 



The map l in the last example has a 'bonus' invariance, namely invariance under transposition: 
= i{A) for any A G GL{3). t is, in fact, invariant with respect to the action of the group 
G = ^2 K (0(3) X 0(3)), with semi-direct product structure 

(r U V)(r' U' V) = I (-'^UU',VV') r = l 

This additional symmetry plays a crucial role in the analysis of the Riemann ellipsoids [381 EJ [12] , 
which are discussed in ^ If we identify so(3) with as usual, and q with x M^, the adjoint 
action satisfies 



Ad(.,f;,V')(r?,C)-| ^u^^vrj) 



(Ut], VC) T = l 

r = -l ' 
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with associated Lie bracket 

[{vX),{fl,C)] = {v^riXxC)- 
A left action of G on SL{3) is determined by 

{1,U,V) ■ A = UAV'^ and (-1, 1, 1) • ^ = 

Unrecognized or unintentional symmetries of invariant maps can lead to unexpected degeneracies 
of the differentials. If such symmetries are found, it is necessary to decide whether they are, in 
fact, relevant to the system or systems at hand; if not, it may be advisable to modify the map 
L SO as to break this extra symmetry. 

4 Relative critical points 

A simplifying feature of the examples treated in y3.1l and l3.2l is the commutivity of the symmetry 
groups. In those examples, the discrete symmetries — negation and reflections across vertical 
planes — determined the nontrivial isotropy subgroups other than the full symmetry group. In 
general, locked Lagrangians, energy-momentum functions, etc. are only invariant with respect 
to the action of the isotropy subgroup of the parametrizing algebra element under the adjoint 
action (or coadjoint action if parametrized by g*). However, the action of the full symmetry 
group typically plays a crucial role in the analysis of the relative critical points. If the (dual) 
algebra elements are treated as fixed parameters, invariant function theory cannot be fully 
utilized. Hence we extend the group action using the (co)adjoint action to obtain a fully invariant 
function on M x g (or M x g*), but only require criticality with respect to the original variables. 

Consider a manifold M with a left G action; assume that the action of the group G on g 
commutes with the adjoint action of G, and that : M x g — >■ R is invariant under the action 

{g, h) ■ (m, ^) = {g-m,h- Adg^) 

of G X G on P = M X Q. (g can be replaced with the dual of the algebra g* and the adjoint 
action replaced with the coadjoint action throughout.) We say that (m,^) is a relative critical 
point of -F iff m is a critical point of the function : M — > M given by F^(m) := F{m$,).. If 
F = f o L for some invariant map t : M x g — > Q and function / : Q — )• M, then (m, ^) is a 
relative critical point of F iff 

d,{m4)f ^'^ev[dl,i^], where L^{m) := i{m,$). (20) 

To illustrate the approach, we begin with a simple and familiar example, the free rigid body. 

Example. 

The free rigid body has Lagrangian L : TS0{3) w 50(3) x M'^ — t- M given by 

L{An) = f{fl) := ifi^m, 

where " : M'^ — > so(3) ~ {3 x 3 skew matrices} is given hy rjy = rj xy for all y € and the 
inertia tensor I is a symmetric positive definite 3x3 matrix. The Lagrangian is invariant with 
respect to left multiplication on 50(3); the associated infinitesimal generator is = 
The locked Lagrangian is thus 

L^{A) = = L{U) = HAA^H) = /(A^O- 
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Since the adjoint action has the form Ad[/^ = when so(3) is identified with M^, the map 
L : SO{3) X ^ given by i{A,^) := is invariant with respect to the action of S0{3) 
on 50(3) X M^, and = f o l. Mechanically, l maps the pair [A, ^) to the body representation 
of ^j^{A). If we identify TaSO{3) with using the left trivialization, then 

di{A,^)iArj,0) = -vA^^ = i{A,^) x rj = [{A^)r^ 

implies that d\i{A,$,) = — with kernel 

kei[d\,[A, ^)] = span{.(A, i)} if ^ / 0. 
Combining this with df{ft) = Ifi yields the relative critical point condition that 

Ii{A,^) G span{/,(^,^)}, 

which is the familiar relative equilibrium condition that the body angular velocity l{A, ^) be an 
eigenvector of I. □ 

The free action and underlying Lie group structure make the rigid body amenable to analysis 
by a wide variety of approaches (See, e.g., [5l[30].) Note that in the treatment given above, there 
is no need to use reduction to obtain the Lie-Poisson structure on M.^ ?a so(3). (Of course, the Lie- 
Poisson structure for the free rigid body is probably more widely known than the corresponding 
unreduced system on TSO{3), but the same strategy can be applied in situations where the 
reduced Poisson is not readily available.) 

The interaction between M and g can be clarified by choosing manifolds Qm, Qgi and Q^, a 
G-adequate map </) : M — > Qj\/ and G x G-adequate map a : g — > Qg, and then seeking a map 
X : M X g ^ such that the map t : M x g — ;> Q := Qm x Q^^ x Qg given by 

i{m,$):={cb{m),x{m,^),am (21) 

is G x G-adequate for M x g. Note that if points (m, ^) with Gm 2 G^ exist, then when 
constructing an adequate map, it is essential to choose x with sufficiently high rank at (m,^) 
to compensate for the possible drop in rank of dm4> relative to nearby points. Given a scalar 
function / on Q = Qm x x Qg, define ^{p,x,a) € T*Qm and ^{p,x,a) £ T*Q^ by 

df df 

— {p,x,a){6p) = df{p,x,a){Sp, 0,0) and —{p,x,a){6x) = df{p,x,a){0,5x,0). 
o<P ox 

If L has the form ()2ip . then (m, ^) is a relative critical point of F = / o iff 

+CX« (|^(^(m,0)) =0, (22) 

where x^i^) •= /t("^)^) fo^ all m € M. Clearly (j22|) is satisfied if 

||(.(m,^)) € ker and ^ ker ■ (23) 

Additional solutions, with ^{i^im,^)) ker [dj^^], exist only if range [d^i?!)] and range [d^x^] 
have nontrivial intersection. 
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4.1 Rigid conditions 

A striking feature of the free rigid body is the coincidence of the configuration manifold and 
symmetry group; the codomain of the invariant map i is the Lie algebra of the symmetry 
group. In general, not all relevant information about a system is captured by the symmetry 
group. However, consideration of the criticality conditions arising from variations tangent to 
the group orbits in M x g can be an important step in the derivation and interpretation of 
the full criticality conditions. The coupling between the given action on M and the coadjoint 
action on g (or g*) induced by the invariance of the function or map often yields significant 
computational simplifications. 

An invariant function F : P satisfies 

dF(m,^)(?7A/(m),ad^^) = 

for any G 0; hence dF^(m)|g.m = iff 



where ^(m, ^) € g* is given by ^{m, $)-r] := dF{m, $){0, rj). (Traditional partial notation will 
be used analogously in what follows.) We will refer to (I24p as the rigid condition. If F = f o l, 
then dM!) is satisfied iff 



Note that changing the dependence of / on (p, or changing the map (p itself, does not change the 
rigid condition. Hence this condition can be particularly useful when studying the relative critical 
points of families of functions with diverse purely M-dependent components. For example, 
changing the potential energy in an invariant natural mechanical system does not alter the 
associated rigid condition. 

Examples. 

• Energy-momentum method 

Consider a symplectic manifold M with G-invariant Hamiltonian H and Ad*-equivariant 
momentum map J : M — )■ g*. Set -F(m,^) = H{m) — j^{m), where : M — )• M is given 
by j^{m) = J{m) ■ ^. The function F is linear in ^, with ^{m,^) = —J{m). Hence the 
rigid condition is ad| J(m) = 0. (See, e.g., [40l Eg [Ml US].) 

• Locked Lagrangian 

Given a G-invariant Lagrangian L on TQ, with a lifted G action, define the locked momen- 
tum map by := J o ^q, where J is the momentum map associated to the lifted action; 
equivalently, 



ad| — (m,^) = 0, 



(24) 




(25) 




(26) 



mq)-C=l-,L^+ec{Q)\ 



6 = ■ 



The rigid condition for F{q,^) = L^{q) = L{^Q{q)) is ad|I^(g) = 0. (See [23].) 
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• Free rigid body 

The invariant map l is linear with respect to the algebra element; hence the locked mo- 
mentum map satisfies 

• C = I f{i{A, ^ + 6 C))|,=o = ^(^' CfHA ^) = C'^AItiA, 

and %(^) = AIl{A,^). It follows that the rigid condition is 

= (AliiA,^)) X ^ = A{IiiA,^)) X LiA,^)). 

Since A is invertible, we recapture the criticality condition that l{A, ^) is an eigenvector 
of I. 

4.2 Internal conditions 

If a pair (m, ^) satisfies the rigid condition for a map F, the relative criticality of the pair can 
be determined by considering the restriction of dF^[m) to a complement Vint to g • m in TmM. 
Criticality conditions complementary to the rigid condition are called internal conditions. A 
convenient choice of complement can facilitate the analysis, but convenience is typically context- 
dependent. In some approaches, the tangent space of an appropriate (local) submanifold can be 
taken as Vint in other situations, Vint is specified via a system of linear equations on TmM. 

If M is a principal G-bundle with projection n : M B, then given a neighborhood U of 
7r(m) and a local section <t : Z// — > V C M, we can take Vnt = d^{m)'^(.'^TT{m)^)'i ("t-)^) is a 
relative critical point of F iff (m, ^) satisfies the rigid conditions and 7r(m) is a critical point of 
the function b i— >■ F{cr{b),(,). If the action is not free, the image of a local section can often be 
replaced by a slice. (See, e.g., [7] and references therein.) A submanifold S of M is a slice at 
m € 5 if 



S is invariant under G 



rri) 



the union G ■ S of all orbits intersecting S is an open neighborhood of the orbit G ■ m; 

there is an equivariant mapping n : G ■ S ^ G ■ m such that 'K\G-m is the identity on G • m 
and ■K~^{m) = S. 

If the rigid condition is satisfied at m and 5 is a slice through (m, ^) then (m, ^) is a relative 
critical point if m is a critical point of the scalar function s i— )■ F[s, $,) on S. 

li F = f o L for some invariant map tand scalar function /, and S is a slice through a 
pair (m, ^) satisfying the rigid conditions, then internal criticality conditions can be obtained 
by replacing the map with the restriction of to S in (j20p . reducing the number of linear 
conditions on df{L{m,^)) from dim M to dim S. Analogous to (j22p and (|23p . (m,^) satisfies the 
internal conditions if 

d*m{cP\s) (|^(^(rn,0)) +<,(xds) (I^W"^'^))) =0- 

and solutions with ^(i(fTi, ^)) ker [dm{4>\s)] exist if 

range [<, (01 5)] n range (^(x^ls)) k^gKO / {0}, (27) 

where /Crig("i,^) ^ -^x(™-^)^^' ^"^^^^"^ is the space of partial derivatives with respect to 

X satisfying the rigid conditions. 
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A rigid-internal decomposition for the energy-momentum method on cotangent bundles of 
principal bundles was introduced in [IQ], and subsequently adapted to the reduced energy- 
momentum method in [H]. An analogous decompositions for locked Lagrangians, including 
systems with nonfree actions, was developed in [23]. All of these decompositions utilize the 
linearization of the rigid condition with respect to m and block diagonalize the second variation 
of the relevant function. Thus these choices of a space of internal variations are well-suited for 
stability and bifurcation analyses, particularly in situations where an appropriate slice is not 
known. 

5 The symmetries of the Riemann elUpsoids 

The Riemann ellipsoids provide a historically significant example with nonabelian symmetry, and 
hence nontrivial rigid component, and diverse assortment of isotropy subgroups. The Riemann 
ellipsoids are steady motions of an incompressible inviscid fluid moving under a self-gravitational 
potential that are generated by one parameter subgroups of the action of SO{3) on the left and 
the group of volume preserving diffeomorphisms on the right [S^. (See [6] for a survey of the 
history of the Riemann ellipsoids and a modern derivation using moments. See |24l [T2] for an 
analysis using geometric methods.) 

Rigidly rotating ellipsoids were used as qualitative models of spinning fluid masses by Newton, 
Maclaurin, and Jacobi [H]. Dirichlet [TT] rigorously proved that affine motions of an ellipsoidal 
fluid mass form an invariant subsystem of the Euler equations for an incompressible inviscid fluid 
under the influence of a self-gravitational potential. This subsystem can be characterized as a 
Lagrangian — or canonical Hamiltonian system on 5*^(3); the right action of volume preserving 
diffeomorphisms on the original configuration manifold determines an action of S0{3) by right 
multiplication. The affine fluid model has configuration manifold SL{3) and symmetry group 
Z2 >< (S0(3) X SO{3)). 

The analysis of the relative equilibria of this system traces the development of symmetry 
methods in conservative systems. Maclaurin analyzed the rigid rotations of axisymmetric fluid 
masses about the axis of symmetry; Jacobi extended these results to triaxial fluid masses [6]. 
Dedekind [10] observed that rigidly rotating steady motions, generated by left multiplication, are 
transformed by transposition to internal swirling steady motions, in which the region occupied 
by the fluid remains unchanged. This transposition invariance played a crucial role in the 
evolution of symmetry-based analyses of dynamical systems, extending techniques developed for 
readily visualizable, intuitively clear spatial rotations to novel group actions — the full body of 
results for the left action immediately yielded a corresponding body of results for the action of 
SO (3) by right multiplication associated to particle relabeling invariance. Riemann [38] built on 
Dedekind's insight, classifying all relative equilibria — including large families of quasi-periodic 
motions — generated by combinations of left and right multiplication. 

We construct an invariant map for the action of G = Z2 k (5*0(3) x 50(3)) on SL{3) by 
restricting the action of Z2 ix (0(3) x 0(3)) on GL(3) described in the last example in §3.31 
and the remarks following that example. Given the relative complexity of this system, we do 
not rigorously investigate the adequacy of this map here, noting only that it satisfies (fTSj) at 
generic points within some key isotropy classes. We consider two classes of points {A, ri, ^) € 
P = SL{3) X X M"^ studied by Riemann. S type points are those for which both components 
of the algebra element are parallel to a common axis determined by A; specifically, an S type 
point is one for which there is a singular value decomposition A = Udiag[a\V^ , U,V ^ 0{3), 
such that 77 = rjUej and C = C^^j foi' some j € {1,2,3} and rj,(^ (^M. Riemann's type I, II, 
and III ellipsoids all have coplanar algebra elements: there is a singular value decomposition 
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A = Udiag[a\V such that U rj and V C belong to a common coordinate plane. As we shall 
see, these states have nontrivial isotropy that reduces the number of criticality conditions to be 
satisfied. 

UAe SL{3) and C := A, then 

ll^ll^ = traceC and ||^~^||^ = ^ ((traceC^ - traceC^) , 

where C := A^ A. Hence 

ct>{A) := {\\A\\W\A~^f) (28) 

can be regarded as the restriction of the invariant map defined in the previous example, with 
the third component — which is identically equal to one on SL{3) — discarded. ([28]) satisfies 

rank[(iA0] + 1 = # of distinct singular values of A. 

The map 

"(^,C):=(hf +||Cf ,(Nf-||Cf)') (29) 

is adequate for the adjoint action on g ~ x M^. (Note that the map (t7,C) ^ll^lP ) IICII^^ 
is not invariant under the adjoint action (jl9p of group elements involving transposition, since 
that action exchanges rj and ^.) 

We now choose a map x depending nontrivially on both M and g. At generic points 
(^{A,r)X) ^ ^) rank four; however, if A has exactly two distinct singular values-and hence 
is rank one-but the isotropy subgroup of (ry, C) has trivial intersection with that of A, then 
X X X a cannot be adequate at {A, rj, ^) unless x has codomain of dimension greater than four. 

The map ()28p is equivariant with respect to the Z2 actions 

A i-^ A~^ and {x,y) {y,x). 

l{A) is fixed iff A posses a unit singular value. In light of this equivariance, and the fact that 
the entries of A~^ are polynomials in the entries of A for A € SL{3), we will use inversion in 
the construction of the invariant map t. Define 

K{A,rjX) :=i||(r/,C)M(A)f = i 

CiA,r,X) := (^7 x AA^rjfAC+{C x ^^^C)' A^rj, 
and 

X{A, r,, C) := {K{A, rj, C),K{A-^, rj, C),V^AC, v^A'^C C{A, rj, C)) . (30) 

The function K and weighted inner products are invariant with respect to negation of the 
algebra element, while C is equivariant with respect to this action and the combination of 
inversion and transposition of A, specifically 

C{A, -V, -C) = -C{A, T), C) = C{A-^, V, C). 

Hence if we define a Z| action on P by 

(-1, 1) • {A, rj, C) = {A, -rj, -C) and (1, -1) • {A, rj, C) = (^"^, V, C), (31) 
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and define a linear Z| action on MJ such that (—1,1) negates the final component, leaving 
the other components unchanged, and (1,-1) pairwise exchanges the first six components and 
negates the seventh component, then the map t = cj) x x ^ Oi determined by (p8]) . ([30]) . and 
()29p is equivariant with respect to this action. This equivariance simplifies the analysis of some 
important classes of relative critical points. 



5.1 Coplanar states 

The time reversal symmetry of the Lagrangian for the affine fluid system leads to invariance 
of the locked Lagrangian under negation of the algebra element. The interaction between this 
action and the G action on P plays an essential role in Riemann's analysis of the affine liquid 
drop system 08]. The invariance of the map l with respect to the G action and equivariance 
with respect to the action ()3ip determine much of the structure of the differential of i at 
points with nontrivial isotropy with respect to these actions. 

We will say that {A, rj, (^) is coplanar if ij and both lie in plane determined by A, i.e. if there 
exist U,V £ SO{3) such that UAV'^ is diagonal, and span{J7?7, VC} is the Euclidean orthogonal 
complement of one of the Euclidean axes. Riemann's types I, II, and III classes of relative 
equilibria are coplanar [38] . The coplanar states are fixed by the two element subgroup of the 
enlarged group G x generated by reflection across the plane containing both components of 
the algebra element. Specifically, given a € with 010203 = 1, ^7, C ^ and U,V € SO{3), 
if we define 

ip{a,T],C,U,V) := {Udiag[a]V^,Ur],VC) and gj{U,V) := {l,U RjU^ ,V RjV^) G G, 

where Rj denotes rotation about ej through vr, then 

gj{U, V) ■ ^{a, rj, C, U, V) = (t){a, rj, C, URj, VRj). 

In particular, if r] and C are both orthogonal to ej, then {gj{U, V), (—1, 1)) fixes ^p{a, 77, C, U, V). 

The equivariance of C with respect to negation of the algebra elements implies that coplanar 
points are both zeroes and critical points of C: if p = ^p{A, r], C, U, V) is coplanar and iIj : P ^ P 
denotes the action of ((7j(?7, y), (— 1, 1)), then C o ip = —C implies —C{p) = C{ip{p)) = C{p) 
and 



dpip*dC{p) = dpip*dC{ip{p)) = d{C o ^/j){p) = -dC{p). 

Since dpip + 1 is nonzero, it follows that p is a critical point of C. 

For notational convenience, assume that p = (/9(diag[a], 77, (^), with ^73 = C3 = Oj the general 
coplanar case follows from invariance. Let (rj,0j/2) be polar coordinates of {rij,(^j), j = 1,2. 
(The motivation for this angle convention will become clear.) We first consider the rigid condi- 
tions: If oi / 02, then the space /Crig(p) of solutions of the rigid condition ([26|) satisfies 

/Crig(p) 5 span{v(a,6'i,6'2), v(a,6li,6'2),e5} , (32) 

where 



v(a, 91,92) := (cos 62 - cos ^1) (iei + 2 64) - w(l, 2) + w(2, 1) 



(33) 



:= cos^i sm9j I aj 63 64 

03 
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and v(a, 9i, 62) is constructed by replacing with — , i = 1, 2, 3, in v(a, 61, 62) and exchanging 
the first and second, and the third and fourth, components of the resulting expression. If 
diag[a] has three distinct singular values, then g • diag[a] is six dimensional and we can take the 
submanifold of diagonal elements of SL{3) as a slice through diag[a]. There exist matrix- valued 
functions N,N : — > M^^^ such that the internal criticality condition is satisfied iff 

= [rl (iV(a) + Nia) smO,) + 1 o 2) (34) 

where by 1 O 2 we mean that the first and second components of a are exchanged, ri is replaced 
by r2, and 9i replaced by 62. Thus in this case, given ^{l{p)) € /Crig(p)> there is a unique value 
of ^{l{p)) compatible with relative criticality of p for f o l. The fifth columns of N{a) and 
N{a) equal zero; hence (0,65) € ker[(i^t(p)] if p is coplanar. 
If ai = 02 = a, and hence 03 = a~^, then 

/Crig(p) 5 spanjei +464,62 + 463,63 - 0^64,65} . (35) 

If we take 



ej := (6j,0,0), i = l,2, and := (0, 6j, 0), z = l,...,5, 

as basis vectors for T*^^Q, then the internal criticality condition is satisfied if df{i{p)) lies in 

df{i{p)) G ker[d^t(p)] = spanja^ei - 62,^(0,77,^)61 + ei + 4e4 + a^(e2 + 4 63), 65} , 

where 

w{a,r]X) :=-(l + a6) (||r,f + ||C|'2^ 



Thus (|27|) is satisfied in this case: range 



n range 



{d;xs) 



l^rig(p) 



is one dimensional 



for generic values of a; hence there are relative critical points with nonzero partial derivative 
with respect to x- Note that in both cases — two and three distinct singular values — ker[(i^i(p)] 
is three dimensional and contains 65, although the spaces of partial derivatives satisfying the 
rigid criticality conditions are of different dimensions. 

Remark. To recapture the rigid conditions for the locked Lagrangian for the actual Riemann 
ellipsoid problem, with algebra-dependent component K{A, t], ^), we find conditions under which 
61 € /Crig(p). Considering p2|) and (p3|) . we see that solutions of the equations 

ai cos 62 sin 61 = 02 cos 61 sin 62 



02 cos 6*2 sin ^1 — ai cos 01 sin 02 = 2 03 (cos 02 — cos 0i) 7^ 



determine states satisfying the rigid condition for the locked Lagrangian. If appropriate non- 
degeneracy conditions on 9i and 62 are satisfied, these equations have the unique solution 



Ol 

03 



2 cos 01 sin 02 02 2 cos 02 sin 0i 



cos 01 + cos 02 



03 



cos 01 + COS 02 



and 



3 (cos 01 + COS 02)^ 

~ sin 201 sin 202 



□ 
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5.2 S type states 

We briefly consider the S type elhpsoids, for which ij and ^ both he on the same Euchdean 
coordinate axis after some change of basis that diagonalizes A. li p = ip{a,7] ejX ^j,U,V) for 
a € A, r],( & 3 & {1; 2, 3}, and U,V G 50(3), then p satisfies the rigid condition, and hence 
is a relative critical point if it satisfies the internal conditions. 
If the absolute values of the entries of a are distinct, then 

ker[d*t] 2 span|vi2 (^-^^ , V2i(a3),sin6'(a3ei - 62) + 40163,63 + 0^64, 65! , 
where (r/, () = (r cos |, r sin |) and Vjj(n) := cie^ + C2(n)ei + cs{u)ej for 

{al - aj) {aj - al) 3 oi as 1 2 ■ n 

ci := 75-^ -, C2[u) := u \-smU, 03(11):= u sm6. 

r^'as 02 ai u 

If p has continuous isotropy component 5^, say p = <f ((a,a,a~^) , 63, C 63, [/, , \a\ 7^ 1, 
then dpi is rank one, with 

ker[(i*t] 2 span{(7? - C)^ei + cei,??^^! + 2ce4,a^ei - 62,61 + 0^62, a'' 63 + 64,65} , 
where c := — 1. 
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